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Ionization-induced  generation  of  strong  Langmuir  waves  by  high-intensity 

Bessel  beam 

V.B.  Gildenburg  and  N.V.  Vvedenskii 

Institute  of  Applied  Physics,  Russian  Academy  of  Sciences,  46  Ul’yanov  Street,  Nizhny  Novgorod,  603950,  Russia 

Dynamics  of  the  field  and  plasma  at  the  gas  breakdown  produced  by 
the  Bessel  (axicon-focused)  laser  beam  is  studied  with  the  effects  of 
natural  plasma  waves  generation  and  reflected  wave  spectrum  con- 
version taken  into  account 


We  study  in  this  paper  the  dynamics  of  the  field  and 
plasma  in  the  discharge  column  created  by  axicon  fo- 
cusing [1-3]  of  a high-intensity  laser  radiation  in  homo- 
geneous gas.  We  focus  on  the  effect  of  natural  plasma 
waves  generation  considered  previously  for  other  types 
of  the  optical  and  microwave  discharges  [4,  5]. 

The  spatiotemporal  evolution  of  the  electric  field 

E = Re[E(r,  f)exP(~*®  01  is  calculated  numerically 
based  on  the  vector  wave  equation  for  the  slow  time 
envelope  of  the  field  E (r,  /) : 

-~+82V(y-E)+zE — V[Vx[VxE]]+rE  = 0(l) 

» 3/  kl 

This  equation  takes  into  account  the  spatial  and  time 
dispersion  and  allows  us  to  describe  the  processes  of 
resonance  excitation  and  Landau  damping  of  Langmuir 
waves  in  the  time-varying  plasma.  In  Eq.  (1),  T is  the 
model  dissipation  operator:  r E = - / a 8 V(V  * E) , e - 
1 - ( N / AC)(I  - rv  / co)  is  the  complex  permittivity  of 

the  plasma,  Nc  = w(co2  + v2)/4ne2  is  the  critical 
density,  V is  the  electron  collision  frequency,  a is  the 
coefficient  of  order  unity,  k$  = co/c , 8 = ^[3  Vp  I co , 
Vp  « c is  the  thermal  electron  velocity. 

Gas  breakdown  is  caused  by  the  optical-field-induced 
ionization  processes.  As  an  example  we  consider  tunnel 
ionization  of  hydrogen  atoms  and  determine  the  time 
averaged  ionization  rate  by  the  known  expression 


Here  Ea  = 5.14xl09  V/cm  and  4.16xl016s'1  are 
the  atomic  field  and  frequency  units,  and  Ng  is  the 

concentration  of  neutral  atoms  before  the  process  of 
ionization.  The  parameter  range  we  are  interested  in 
here  is  the  following:  the  wavelength  A ~ 0.6-10  |im, 

laser  pulse  intensity  S ~ 1014 -1016  W/cm2,  and  the 
gas  pressure  p ~ 0.3  - 70  atm. 

We  consider  the  model  of  the  axially  symmetric  dis- 
charge: N (r,  t)  = N ( r , /)  produced  by  the  rotating  cy- 
lindrical wave  with  the  complex  envelope  of  electric 
field  E(r,  t)  = E (r,  <p,  z,  t)  = E (r,  t)  exp(/cp  + ik$  z cos  0) 
(r,q>,z  are  the  cylindrical  coordinates).  Outside  the 


plasma  (r>  R,  N(r  > R)  = 0 ) the  field  is  a superposi- 
tion of-  the  converging  (incident)  and  diverging  (re- 
flected) TE  and  TM  waves  with  a given  angle  of  incli- 
nation 0 to  the  axis  of  symmetry  z . The  incident  wave 
is  given  by  the  axial  components  of  electric  and  mag- 

netic  fields:  E{Jn)  =C(t)Hx(2)  (k0  r sin  Q)  y H(Jn)  = 

-icosQE{Jn),  C(t)  = Aexp(-(t-t0)2/x2),  //,(2)  is 

the  first-order  Hankel  function  describing  the  converg- 
ing wave.  The  correlation  between  the  amplitudes  of 
these  components  (the  coefficient  — / cos  0 ) is  chosen  so 
that  the  transversal  components  of  the  fields  in  the  ab- 
sence of  plasma  are  circular  polarized  and  are  the  zero- 
order  Bessel  function  of  radius  r : E^(r)  -iEr  (r)~ 

J o (k$  r sin  0).  We  note  that,  in  the  presence  of  plasma, 
the  field  remains  circularly  polarized  only  on  the  axis. 

Equations  (1)  and  (2)  were  solved  numerically  in  the 
space  interval  0 < r < R with  the  initial  conditions: 
N(r  ,0)  = 0 , Ey  (r,  0)  = -2  C (0)  cot  0 J0  (^or  s*n  = 

iEr (r,0) , Ez  (/%0)  = 2C(0)</1  (#o  rsinO)  and  the 

following  boundary  conditions:  (i)  the  solution  is  ana- 
lytical at  r = 0:  Ez  = 0,  d Er  ^ ! dr= 0,  (ii)  all  field 

components  (including  Er  ) are  continuos  at  r - R (for 

detailed  expressions  and  explanation  see  [6]). 

It  has  been  found  that  the  scenario  of  the  breakdown 

process  depends  greatly  on  the  convergence  angle  of  the 

wave.  If  this  angle  is  less  than  some  critical  value 

0 2 
0C  ~ 25  , the  maximum  plasma  density  Nmax  ~ Ac0 

that  is  less  than  the  critical  one.  However,  at  the  angle 
exceeding  the  critical  one,  the  plasma  density  at  the  axis 
increases  in  the  sharpening  regime  and  passes  the  criti- 
cal point,  after  that  the  fast  ionization  wave  containing 
the  plasma  resonance  point  at  the  leading  front  propa- 
gates in  the  radial  direction.  The  results  of  numerical 
calculations  are  presented  on  Figs.  1 and  2 in  dimen- 
sionless variables  k$r  -» r , cot-W,  E/£fl-*E, 

N/Nc=n  for  the  parameter  values:  0 = 30°, 

£2/(0=  22,  *o8=  y[^vT /c  = 0.02,  v/(o  = 0.01, 

a = 0.1,  Ng-l.5Nc , /0G)=100,  tco=50, 

Al  Ea=  0.0204  , k o R = 4 . At  the  given  values  of  A , 
the  maximum  field  at  the  axis  in  the  absence  of  plasma 
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Fig.  I.  Evolution  of  the  plasma  density  n(r,t ) and  the  field 
amplitude  | E (r,/)  | at  0 = 30° . Curves  1-9  correspond  to  the 
time  instants  /=  80,  90,  95,  105,  110,  120,  130,  150,  and 
200,  respectively.  The  dashed  curve  shows  the  unperturbed 
field  | Evgc(r)  | (in  the  absence  of  a plasma)  at  / * t0  = 100. 

is  the  same:  | E |max  / Ea  = 0. 1 . The  above  dimension- 
less parameters  correspond  to  the  vacuum  wavelength 
^.“0.8  |Im,  maximum  pulse  intensity 

5,«3xl014  W/cm2,  pulse  duration  (at  the  level  of 
Me)  t^~2~  ~ 30  fs,  and  the  gas  pressure  p ~ 60  atm . 


Fig.  2.  Time  behavior  of  the  amplitude  function  R cG{t)  of 
the  reflected  TE  wave  at  large  times.  The  curve  illustrates  the 
radiation  from  the  discharge  at  the  upshifted  frequency 

^pmz\  ~ 1 *22(0  . 


the  scattered  spectrum  of  an  ionizing  electromagnetic 
wave  (previously  described  in  the  model  of  a thin  gas 
slab  [5])  is  illustrated  in  Fig.  2 by  the  time  dependence 
of  the  quantity  ReG(/)  determining  the  amplitude  of 

H~  -component  of  the  reflected  wave  = 

{kyr)  at  large  times  t > 220  (after  the  end  of 

the  incident  pulse),  when  the  signal  amplitude  at  the 
fundamental  frequency  has  already  substantially  de- 
creased, but  the  Langmuir  oscillations  still  exist  and 
continue  emitting.  The  transversal  wave  number  at  the 
frequency  0^^  is  *±  = («pmax  /c)sin  9p  = 

^((flpma x/c)2  -^o  c°s2  9 > that  is  the  corresponding 

reflected  cylindrical  wave  is  inclined  to  z-axis  at  the 
angle  0^>0  (cosO^  = (to  / to  p max)  cos  0 ).  In  the 


The  transition  of  the  plasma  density  through  the  criti- 
cal value  at  0>  0C  (Fig.  1)  is  accompanied  by  the  exci- 
tation of  intense  Langmuir  oscillations,  whose  ampli- 
tude reaches  its  maximum  (twice  as  high  as  the  ampli- 
tude of  the  unperturbed  electric  field  at  the  axis)  at  the 
front  of  the  ionization  wave  at  r = 1 . 

The  coupling  of  the  excited  Langmuir  oscillations  to 
an  external  electromagnetic  field  (due  to  the  presence  of 
a fairly  sharp  boundary  of  the  ionized  region)  gives  rise 
to  the  partial  emission  of  their  energy  into  the  sur- 
rounding space,  i.e.,  to  the  occurrence  (along  with  the 
fundamental  frequency  component  w ) of  one  or  several 
components  at  frequencies  close  to  G)pmax  = 

4iie2N g / m ~ 1.22o)  in  the  spectrum  of  the  cylindri- 
cal wave  reflected  from  the  discharge.  In  view  of  the 
linear  character  of  the  "transition"  resonant  excitation  of 
Langmuir  oscillations,  the  intensities  of  the  shifted 
spectral  components  are  proportional  to  the  intensity  of 
the  incident  wave.  This  linear  parametric  conversion  of 


given  numerical  example,  0 = 30°,  0^=45°;  at  t~ 
250  the  intensity  of  the  frequency  shifted  component  is 

about  lO”4  of  the  maximum  intensity  of  the  incident 
wave. 
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